Chaos in the thermal regime for pinned manifolds via functional RG 
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The statistical correlations of two copies of a d- dimensional elastic manifold embedded in slightly 
different frozen disorder are studied using the Functional Renormalization Group to one-loop, order 
0(e = 4 — d), accuracy. Determining the initial (short scale) growth of mutual correlations, i.e. chaos 
exponents, requires control of a system of coupled differential (FRG) equations (for the renormalized 
mutual and self disorder correlators) in a very delicate boundary layer regime. Some progress is 
achieved at non-zero temperature T > 0, where linear analysis can be used. A growth exponent a 
is defined from center of mass fluctuations in a quadratic potential. In the case where temperature 
is marginal, e.g. a periodic manifold in d = 2, we demonstrate analytically and numerically that 
a = e(l/3 — 1/(2 ln(l/T)) with interesting and unexpected logarithmic corrections at low T. For 
short range (random bond) disorder our analysis indicates that a — 0.083346(6)e with large finite 
size corrections. 



I. INTRODUCTION 



A. Overview 



Systems with quenched disorder are especially sensi- 
tive to small external perturbations. This phenomenon- 
is called chaos. Upon a perturbation of amplitude S (an 
energy scale), e.g. a small change in temperature (tem- 
perature chaos) or in the disorder (disorder chaos) the 
configuration of the system, while only weakly affected 
at small scales, changes completely beyond a length scale 
L$. This overlap length diverges as Lg ~ 5~ a for small 6, 
a being called the chaos exponent. Chaos can be studied 
either at T = via the sensitivity of the ground state to 
perturbations, or in the thermal regime at T > 0. Chaos 
has been studied in spin glasses and other disordered sys- 
tems using droplet arguments^* 3 -, numerics^ and mean 
field calculations^ . While central to disordered systems, 
chaos is still not fully understood. There has been some 
controversy as to whether the overlap length is finite or 
infinite (no chaos), in high dimensions as well as in mean 
fieldM. 

An interesting class of systems exhibiting chaos are 
elastic systems in random potentials such as domain 
walls in disordered magnets 7 , or periodic systems such as 
charge density waves^ and vortex lattices on disordered 
substrates*^. They have energy-dominated glass fixed 
points at which temperature is irrelevant, Tl ~ L — , 
6 = d — 2 + 2( being the free energy fluctuation ex- 
ponent. The roughness exponent £ controls the scal- 
ing of the typical deformation u ~ IS with the internal 
size L of the pinned configurations of the elastic object. 
Chaos in pinned manifolds was studied mostly via scaling 
arguments**. The directed polymer (d = 1) was stud- 
ied numerically and via analytical arguments for N = 1 
indicating a = 1/6 in agreement with droplets* 2 -, and 
recently, on hierarchical lattices* 3 -. In d = 2 chaos was 



demonstrated for periodic systems near the glass transi- 
tion T g using the Cardy Ostlund RG— . 

A successful approach to disordered elastic systems is 
the functional renormalization group (FRG). It allows 
an efficient determination of the roughness of the ground 
state within an expansion in e = 4 — d, to one loop 15 ' 16 i 17 
and recently to higher orders**V FRG involves a coupling 
constant function, A(u) which measures the renormalized 
correlator of the pinning force and becomes non-analytic 
beyond the (Larkin) length scale L c where pinning pro- 
duces metastability. FRG was also studied at non-zero 
temperature to one loop**. There it was shown*£ that 
rare thermal fluctuations (droplets) lead to the cusp of 
A(u) being rounded off inside a region of size u ~ T, 
the so-called thermal boundary layer (TBL). Recently 22 
the full function A (it) was shown to be a proper physical 
observable, which describes the fluctuations of the center 
of mass of an elastic manifold confined by an harmonic 
well. It was determined by a high precision numerical 
calculation 23 at T = and found to compare remark- 
ably well with analytical predictions, already at the one 
loop level. The cusp in A(u) was also observed^ 3 -, and 
shown to result from so called shocks, abrupt switches of 
the manifold from one ground state to another, as the 
position of the center of the well is varied. 

When applying FRG to the problem of chaos^, one 
must follow not only the flow of A(u), but also the flow of 
a second FRG correlation function D{u) which encodes 
the mutual correlations between the centers of mass of 
two manifold copies seeing slightly different disorders. 
These flows have been analyzed mostly for large scale 
mutual correlations. It was found2i that the residual 
correlations decay to zero at large scale for (i) the 'ran- 
dom periodic' class (RP), i.e. a correlator A(it) peri- 
odic in u (which describes charge density waves and vor- 
tex lattices) and for (ii) the 'random bond' class (RB), 
i.e. a correlator A(it) = —R"(u) where R(u) is a short 
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range function (which describes magnetic domain walls 
in short range disorder). Hence there is chaos with a 
finite L$. For the 'random field' class (RF), which de- 
scribes magnetic interfaces in the presence of random 
fields, it was found that residual correlations remain non- 
zero. However, determining the initial short scale growth 
of mutual correlations, and hence the chaos exponents, 
is difficult. It requires good control of the system of cou- 
pled differential one-loop FRG equations, specifically of 
the separation — initially very small — between D(u) and 
A (it). Since D(u) remains analytic (which was confirmed 
numerically 23 ) while A(u) develops a cusp at u = 0, one 
cannot use standard linear analysis. Instead, the chaos 
boundary layer, where the two functions differ, must be 
investigated: a non-trivial task in full generality. 

The aim of the present paper is more modest. We 
study the chaos problem on the one-loop FRG equation 
in the thermal regime T > 0, where linear analysis is ap- 
plicable. The easiest case to handle is when temperature 
does not flow under RG and there is a line of fixed points 
as temperature is varied. This happens for the random 
periodic class and d = 2, the case on which we focus here. 
Hence for each T > there is an analytic FP, At (it), and 
one can use linear analysis to extract the growth expo- 
nent for mutual correlations. It is still a non-trivial task 
as one must perform the analysis both inside and out- 
side the TBL, and match the two results. Surprisingly, 
as At(u) becomes non-analytic for T — > 0, one finds 
logarithmic corrections to the growth exponent. These 
corrections are confirmed by a careful numerical study 
of the differential equation. Despite being a special case 
of the full chaos problem, the random periodic class for 
al = 2 already illustrates the difficulty of obtaining the 
accurate behavior. An extension to the random bond 
class is then proposed, again within the thermal regime. 



B. The model and the observables 

We investigate how two identical copies of an harmonic 
elastic manifold embedded in frozen disorder decorrelate 
when they are exposed to slightly different disorder. Here 
we focus on interfaces, i.e. manifold whose deformations 
are parameterized by a real valued displacement field 
u(x), where x is the d-dimensional internal coordinate. 
The system is described by the following Hamiltonian 



Hv,M = j ^ 



(1) 



+ V i {u i (x) 1 x) + -n^iu^x) - v) 2 



where the two copies i = 1, 2 are not mutually interact- 
ing. They are however coupled via the correlations be- 
tween the two random potentials, and between the corre- 
sponding random pinning forces Fi(u,x) — —d u Vi(u,x), 



whose correlation matrices take the form: 

Vi(u, x)Vj(u', x>) = R\f(u - u')5 d (x - x') (2) 
Fi(u,x)Fj(u',x') = Ag 5 (u - u')5 d {x - x 1 ). (3) 



with Ag } ( M ) 



-Rf>"{u). 



We use the superscript 
(0) to denote bare disorder, to distinguish it from the 
renormalized disorder defined below. The form of these 
correlation functions differentiates the three main uni- 
versality classes: (i) random periodic (RP): a periodic 
RS Q >{u) (ii) random bond (RB) a short range function 
R(°>(u) (iii) random field (RF) RW(u) ~ a\u\ and 
A<°)(u) short range. One way of realizing eq. @ is 
to consider two disorder copies of the following form: 
Vi(u,x) = V(u,x) ± 5W(u, x) with V and W mutually 
uncorrelated Gaussian disorders of zero mean, and 5 a 
very small parameter. In the simplest case V and W 
are in the same class with identical correlator. Note that 
a quadratic well, i.e. a mass term, has been added to 
confine each elastic manifold copy to a mean position v, 
another external parameter. Its role is discussed below. 
One is usually interested in the limit m — > 0. 

The statistical correlations of the two ground states 
are measured by the standard observables, i.e. the cor- 
relation functions: 



dj{x - x') = ((u^x) - u i {x'))){{v?[x) - v?(x'))) (4) 

where — denotes disorder averages. We will also study 
non-zero, though low, temperature, denoting the thermal 
averages by (..) (using the canonical partition function 
Zv.v — f Due~ Hv ^ u ^/ T ). For a single copy one expects 
mean square deformations due to disorder to result in 
Cn(x) ~ en |x| 2C at large L c < x < ^, ( being the 
roughness exponent, independent of temperature at low 
T (and in the whole glass phase where the manifold is 
pinned near its ground state with only a few active ther- 
mal excitations). Standard arguments that assume the 
existence of a single diverging scale, the overlap length 
Lg ~ S~ l / a , suggest the following scaling form for the 
two-point correlation function between different copies 
at large L c < x < -|: 



C 12 (x) =x 2 ^(Sx a ). 



(5) 



The overlap length separates correlated from uncorre- 
lated scales, and depends sensitively on the small dif- 
ference 5 between the bare disorders of the two copies. 
In analogy with chaotic dynamical systems, in which tiny 
differences in initial conditions are amplified via the Lya- 
punov exponent(s) to large scale differences, one intro- 
duces the chaos exponent a. This exponent is a measure 
of how the two copies effectively split as scale increases. 
Qualitatively, this splitting is characterized by a dimen- 
sionless scale-dependent parameter which grows under 
RG as Si = Se al where I is the log-scale, e.g. / = InL. 
This behavior is suggested by droplet arguments, first de- 
veloped for spin glasses where they predict^ a = df/2 — 6 
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with —9 being the thermal eigenvalue and df the fractal 
dimension of the droplets. In the case of manifolds the 
same formula was proposed^ii for the SR disorder class, 
with df — d, namely a = d/2 - 8 = (e - 4£)/2. 

Another observable, introduced recently^ 2 -, quantifies 
the fluctuations of the center of mass confined by an har- 
monic well. One defines u l (x;v) — (u l (x)} the thermally 
averaged position. It depends on the position of the cen- 
ter of the harmonic well v. One denotes the center of 
mass of the manifold by u l (v) = L~ d J d d xu l (x;v), L d 
being the system volume. The second cumulant of its 
position, as the disorder is varied: 

rnti&iv) - v)(ui(v') - v') = L _d Ay (v - v') (6) 

defines the renormalized pinning force cross-correlator 
Ajj(u). At zero temperature these functions measure the 
correlations between the shocks in the two copies. These 
abrupt jumps of the manifold as the center v is varied 
do not occur exactly at the same place in the two copies, 
which results in the cross correlator Ai2(u) remaining a 
smooth function of u. This was confirmed in— where si- 
multaneous shocks in the two copies were examined and 
A 12 (it) was computed numerically. Another useful quan- 
tity is the free energy of each copy Vi(v) = —ThxZv,v' 
for each disorder configuration it is a random function of 
the well center position v, hence one defines the renor- 
malized second cumulant of the potential: 



V i (v)V j (v') = L d R ij (v-v') (7) 
and it is easy to see that Ay(ii) = —R"Au). 



C. Functional RG approach 

The functional RG method allows us to compute from 
first principles the observables defined above, namely the 
correlation functions of eq. (j4]) and the renormalized cor- 
relators of eq. © and 0. FRG is based on the repli- 
cated field theory, and proceeds via a loop expansion 
for which — at zero temperature — the small param- 
eter is e = 4 — d. Here we give only results, and refer 
to 22 ' 24 for reviews on the method. The FRG flow of the 
renormalized correlators Rij were derived ir>2i. The FRG 
flow equations for the correlators Rij and Ay , defined in 
eq. ((6]) and 0, are found by computing the effective 
action and its flow equation, as the mass is varied. One 
defines m = mi — moe~ l where I is the usual RG logarith- 
mic scaled. We define the rescaled dimensionless force 
correlator A,j via Ay(u) = A' d 1 m' i [ 2 ^Ajj(itmf), with 
1 = 87r 2 , and the roughness exponent £ reflecting the 
sclf-affine scaling property of the manifold. A rescaled 
temperature is also defined as 7} = Tm\. We simplify 
notation by writing A = An for the one-copy correlator 
and D = A12 for the two-copy correlator. Their flow 



equations are^i, respectively, to one loop: 

diA(u) = (e-2C)A(w) + C«A»- A'(u) 2 (8a) 
+ [A(0)- A(u)]A"(u) + TjA"(u) 

diD{u) = (e-2()D(u) + (uD'(u)-D'(u) 2 (8b) 
+ [A{0) -D(u)]D"(u)+T l D"{u) 

The zero temperature FRG equation are obtained by set- 
ting Ti = 0. 

Note that the temperature is irrelevant (for 9 > 0) 
because Tj = Ter el (setting too = 1). However, the 
TjA"(u) term keeps the correlation functions smooth for 
any non-zero T. If one studies the FRG to only one loop, 
one can also use the (more qualitative) Wilson RG proce- 
dure, which consists of varying the short scale momentum 
cutoff A/ = Aoe~'. In that case the mass cutoff is unnec- 
essary: one can set to = and estimate the correlation 
functions of eq. ((4]) at non-zero momentum. 

It is important to note that the one-copy correlator in 
eq. (|5|) evolves independently, whereas the two-copy cor- 
relator is linked to the former via its value at the origin 
A(0). This small but crucial difference entails opposing 
behaviors for the two correlators^: the one-copy correla- 
tor converges to its stable fix point, whereas the two-copy 
correlator diverges towards another fixed point. The dif- 
ference between the two correlators is denoted 

0, = A(0) - D(0) (9) 
6 i=0 5 d {x-x') = 

7j(-F\(0, x) — -F 2 (0, x))(Fi(0, x 1 ) — F 2 (0, x')) 

One easily sees from eq. (j8]) that 6; generates an addi- 
tional term 0;D"(u) in the flow of D(u), compared to 
that of A (it). Hence one can think of 0/ as an effec- 
tive temperature (its bare value before renormalization 
is ©;=o ~ S 2 ) but its flow under RG is very different 
from that of the real temperature T . In fact, its flow is 
determined self-consistently by the two equations. While 
in general it smoothes the form of D(u) which hence re- 
mains analytic— , it is often relevant, i.e. grows with I, 
by contrast to temperature. 

To compute the observables of eq. ([6]) and as a func- 
tion of to one must solve the above flow equations. The 
difference (eq. ©) is a direct measure of the fluctuations 
of the distance between the (thermally averaged) centers 
of mass u l (v) of the two copies: 

(u l {v) -u 2 {v)) 2 = A d L~ d m- d - 2C Qi (10) 

The correlation function of eq. ([4|) is more delicate to 
compute. The general formula Cij(q) = Ay i ;(0)/(g 2 + 
to 2 ) 2 is exact for q — 0, to = mi, and holds to O(e) (i.e. 
to one loop accuracy) for q ~ to. It also holds for q = A;, 
to = within the (one loop) Wilson scheme and hence 
provides an estimate for the correlation function at large 
scale (and for q^> m): 

C l3 (q) « Aaq-^A^e^ 1 ]^^ (11) 
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However caution is required when using this estimate, 
even when computing the simpler, small q behavior. For 
instance, for the random periodic (RP) universality class 
(for which £ = 0) the estimate is correct only for d > 2. 
In general one must examine more carefully the FRG 
for the non-local part of the effective action 2 ^: in d = 
2, 9 = the two gradient term becomes dominant and 
yields an extra ln(l/<7) in the single copy correlation and 
the famous In 2 |x| in the real space correlation^. The 
study needed to elucidate the initial growth regime being 
even more subtle, our results here will mostly concern the 
q = 0, center-of-mass behavior of eq. (flT))) and ©. 
Finally note that the thermal correlations: 

C%(x-x') = (12) 
JJvMx) - u l {x')){u3{x) - w>(x'))) 

- {{U l (x) - U l (x'))){{ul(x) - Ul(x'))) 

always vanish identically for i j, because the two copies 
do not interact. 

The outline of the paper is as follows. We start by 
studying the simplest case of the random periodic class. 
In section|n]we focus on zero temperature, review known 
results and explain why the problem is difficult. In sec- 
tion IIIII we study the RP problem at T > in the so 
called marginal case of 6 = (i.e. in d = 2) where tem- 
perature does not flow. This allows to use linear analysis. 
In section llVl analytical predictions and numerical analy- 
sis are compared. In section [V] we generalize these results 
to the random bond class, and RP for d = 4 — e. The 
results are summarized and discussed in the conclusion. 



II. RANDOM PERIODIC CLASS: ZERO 
TEMPERATURE CONSIDERATIONS 

A. Flow of single copy and fixed point of the 
random periodic problem (CDW, Bragg glass) 

Here we study the random periodic class (RP), which 
has logarithmic roughness, i.e £ = 0. The single-copy 
correlator is a periodic function A(u + 1) = A(u) of nor- 
malized period one, and obeys the following flow equation 
at T = 0, from eq. ©: 

d t A(u) = eA(w)-A» 2 (13) 
+ (A(0) - A(u))A"(«) 

where e = 4 — d. This flow is well understoo d 15 ! 17 ' 19 . 
Beyond the Larkin scale — here m = m c = moe 1 " - 
it develops a non-analyticity (cusp) at u — 0, and flows 
towards an attractive fixed point: 

A» = ^(l-6u(l-t*)) , u€[0,l] (14) 

Note that this correlator presents a cusp at u = because 
of the periodicity condition A(u + 1) = A(u). 



B. Flow of the two-copy correlator 

The flow of the two-copy correlator is more intricate, 
due to the (scale dependent) coupling 6; to the single- 
copy correlator. 

diD(u) = eD(u)~D'{u) 2 (15) 
+ {Qi+D(Q)-D{u))D"{u) 
= eD{u)-D'(u) 2 

+ (A(0) - D{u))D"(u) 

The only difference between the two equations is the 
term ®iD"(u). If one starts with a very small Go ~ S 2 <C 
1, the two correlators A(w) and D(u) remain practically 
identical up to the Larkin scale l c , and the difference 
0; = A(0) — D(0) remains small. One finds that it grows 
as Qi ~ e 2ajji Oo, with = (e — 2£)/2, since in most 
of this regime one can neglect the non-linearities. Near 
the Larkin scale, non-linearities become important and 
A" becomes large as the cusp develops. Once 8; grows 
such that OiD"(0) - O;A"(0) - eD(0), which occurs 
very near Z c , the two-copy correlator D(u) starts to differ 
from the one-copy correlator. By analogy with temper- 
ature, one expects this difference to be mostly confined 
to a boundary layer (BL) of width u ~ 8; C 1 around 
u = 0. This BL is called chaos BL to distinguish it 
from the thermal BL u ~ T; (absent at zero tempera- 
ture). While the one-copy correlator flows towards its 
fixed point and develops a cusp (eq. (fl4)) ). the two-copy 
correlator remains analytic. 

The flow beyond the Larkin scale is non-trivial. In 
Ref^ the final behavior of the flow for / — > oo was 
examined. It was found that ultimately D{u) flows to 
D(u) — for the RP class, hence there are no residual 
correlations between the two copies. Here we address a 
different question. We are interested in the first phase of 
the FRG flow, i.e. we study how the difference 8; grows 
beyond the Larkin scale. 

Clearly A* is a fixed point both for D(u) and A(u). 
However while A flows towards its attractive fixed point 
A*, D is repelled by it. Let us assume that the 
one-copy correlator has already reached its fixed point 
A(u) = A* (it) (eq. (HD). We then have to solve the 
FRG equations for the flow of the two-copy correlator D, 
from which we can deduce the behavior of the difference 
8 ; = A*(0)-£>(0) : 

diD{u) = eD(u)-D'(u) 2 

+ (A*(0) - D(u))D"(u) (16) 
die t = -a ; ,D(0) 

One can always write: 

£)(«) = A*(u) + /(u,l) (17) 

where, during the initial growth phase, f(u,l) remains 
small (in a sense to be made precise below). The prob- 
lem we face at T — 0, is that the one-copy fixed point 
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correlator is non-analytic, whereas the two-copy corre- 
lator is analytic, with the cusp rounded off inside the 
chaos boundary layer. Hence the function f(u,l) should 
be equally non-analytic, to cancel the fixed point cusp 
and leave a smooth analytic function D(u). 

An important property of f(u, I) arises from the poten- 
tiality constraint. Ay(n) = —R'^(u) can at most have a 
cusp singularity— and for the RP class the Rtj are peri- 
odic, while for the RB class they must be short ranged, 
which implies: 



where 2 Fi(a,/3, 7> z) = E^o i£ feyf k lf> with («)» = 
r(a + n)/r(a), being Gauss' hypergeometric series. It 
provides a convergent series for < u < 1, whose behav- 
ior near u = is: 



fa(u) = A a 



1 



+ B a In«(l + bi(a)u + b 2 {a)u 2 + ..) 



+ Co (a) + ci(a)u + c 2 (a)w 2 + 
with B a = 12- - 6 and 



(24) 



duf(u, = 



(18) 



for both the RP and the RB class (the integration do- 
mains being u £ [0, 1] and u £ [0, oo], respectively.) 

We now attempt a linear expansion around the fixed 
point. We need eigenfunctions f(u,l) that are non- 
analytic and that obey the zero-mean constraint of 
eq. (fT5)) . Although linear analysis is not necessarily valid 
inside the boundary layer, outside (for u 3> O/) it is ap- 
propriate. 



C. Eigenvalue problem of the T — linearized flow 
equation 

One starts from: 

A*(0)-£>(u) = A*(0)-A*(u) + /(u,i) (19) 

= i u (l-u)+f(u,l) (20) 
6 

and inserts it into the flow equation (|16p keeping linear 
terms only. We assume the eigenvector flows as: 



f(u,l) = exp(2aZ)/ a (u) 



(21) 



which provides one definition of a growth exponent a, as 
discussed below. One gets: 

= (e-2o)/(u) + |[«(l -«)/(«)]" (22) 

= 4(l-^)/(«) 

+ 2(1 - 2u)f'(u) + u(l - u)f"(u) 

For any value of a this linear second order differen- 
tial equation has two types of solutions on the interval 
u € [0,1], even and odd about u = 1/2. We must se- 
lect the even one (since the correlators are symmetric: 
A(— u) = A(it) and D(— u) = D{u) combined with peri- 
odicity). The even solution for general eigenvalue a can 
be expressed in terms of the hypergeometric function 2-F1 



./;,(„= ,F\[ 1-1/25-482, 



(23) 



4T[f - 1^/25 - 48f]r[f + \^25~- 



(25) 



Hence for generic a there is a non-integrable divergence 
at u — together with terms non- analytic in Note 
a particularly simple solution for a = e/2: f a ( u ) — 
1/(m(1 — u)). It is however not integrable at u = 0. 

For f a (u) to be defined (i.e. finite) one needs A a = 
0, which gives: 



-± V25-48- = -n 
4 4V e 



(26) 



where n is a positive integer. It yields a series of values: 

71 = 0,1,2,.. (27) 



a = a n 



1 1 2 1 

n n l e 

3 3 2 



for which the hypergeometric series becomes a polyno- 
mial of finite order. For the highest value: 



fa(u) = c 



(28) 



the eigenfunction is a constant. The next one is a = —% 
and corresponds to f a [u) = c e (l — 5it(l — u)). A more 
detailed analysis is performed in appendix [AJ It is found 
that none of these eigenfunctions satisfy the potentiality 
condition of eq. $TE\i f a (u) — 0. 

This dilemma of non-zero mean eigenfunctions will be 
overcome in the next section, by studying the FRG flow 
at non-zero temperature. Setting T > induces the cusp 
to round off within a boundary layer around u = 0, and 
permits solutions of the eigenvalue problem that have 
zero mean value. 

More generally, the above analysis is valid outside the 
boundary layer (BL), be it a thermal BL u ~ T; or the 
chaos BL u ~ G;. One expects the blow-up of the eigen- 
function near the origin to be rounded off within the BL. 
We now examine how. 



III. RANDOM PERIODIC UNIVERSALITY 
CLASS: NON-ZERO TEMPERATURE T > IN 
THE MARGINAL CASE (6 = 0, d = 2). 

To escape from the difficulty of a non-analytic fixed 
point A* we now consider the problem at non-zero tem- 
perature. In this section we focus on the simplest case, 
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6 = d — 2 = 0, where temperature does not flow under 
RG. Hence there is a line of analytic fixed points A^, in- 
dexed by T > 0, which converge to A* = A^ =0 at the end 
of the line T — > + . Around each of these fixed points 
linear analysis is then possible for all u. The physical 
temperature T introduces a cutoff scale for u that allows 
us to find a coherent solution to the eigenvalue prob- 
lem. Of course here e = 2, hence the one- loop results 
are expected to be approximate. Several recent works 
have found that the one-loop scheme provides reasonable 
approximations of exponents and a clear, qualitatively 
correct picture for this mode l 23 ' 26 ' 27 . 



A. One-copy correlator at non-zero temperature T 

Temperature enters the FRG flow equation for the one- 
copy correlator in a natural way: 

diA T {u) = eA T {u) - A' T {uf 

+ [T + A T (Q)-A T (u)}A'+(u) 

The resulting fixed point equation is integrable: 

eA* T {u) = A* T '(u) 2 



(29) 



(30) 



(T + A* T (0)-A* T {u))A* T "{u) 



1, 



= ^[(A^)-A^(O)-T) 2 ]" 



and is implicitly solved by quadraturi 



with 



-G(T,T + A* T (0)-A* T (u)) (31) 



dy y 



V(y - T )(v - y-)(v+ - v) 

3AJ(0) + T 



(32) 
(33) 



± 



3(3A£(0) + T)(A£(0) + 3T) 
with y_ < < T < y + . 



G(a,b) = 
4y± = 

The constraint 



:G(T, y + ) yields the value of the 

fixed point correlator at zero A^(0) as a function of T. 

The finite temperature correlator fixed point A^(u), 
given implicitly by equation (|3 1 1) reduces to the non- 
analytic zero-temperature correlator of eq. (fT4")l as T — > 0. 
Notice how the finite temperature T rounds off the cusp 
within a boundary layer of width ~ T: The curvature at 

i *,,(n\_ an d within 



the origin becomes finite A^"(0) = 
the boundary layer, for u <C T, the following series ex- 
pansion holds: 



A* T {u) = A* T (Q) 



1 - 



A^(0)T 
2 

3A^(0) + 
12 

15A£(0) + T 
30 



(34) 



T |ti\ 2 



T 



O 



(¥)' 



These are the first terms of a systematic low temperature 

1 Q 

expansion : 



(35) 



k>l 



with 4> k (0) = valid inside the TBL, i.e. for u/T = 0(1). 
The following matching conditions hold at large x — u/T: 

M*) ~c k \x\ k (36) 
E fc >x c k \u\ k = A* T=0 (0) - A* T=0 (u) (37) 

since for u = 0(1), A^(u) = A* T=0 (u) +0{T). The first 
scaling form is </>i(x) = <f>(%) — 1 where we define: 



1 



6 



(38) 



as can be verified by inserting the expansion of eq. (|35|) 
into the FRG flow eq. (|2"9"]) and collecting orders in 0(T). 
The low temperature expansion will be detailed and gen- 
eralized in section [Vj where (j>2 will also be computed. 
For the present purpose the following low temperature 
form of the BL is sufficient: 



(39) 




l+(^) 2 -l +T 2 Uu/T) + 0(T*) 



It reproduces well the T — * limit si, as we n as the 
first term of the power series expansion around u = of 
eq. (|34]). 



B. Linearization of non-zero T flow equation 

The FRG flow equation at T > for the two-copy 
correlator, assuming the one-copy correlator has reached 
its fixed point Ay (it), reads: 

diD{u) = eD{u)-D'{u) 2 (40) 
+(T + A* T (0)-D(u))D"(u) 

with di®i = —diD(0). Since D(u) = A^(u) is now a 
(analytic) fixed point of this equation, we define: 



D(u) = A*(u) + f(u,l) 



(41) 



and perform a linear analysis for small /, i.e. we write 



„2aZ 



f a {u) and look for eigenfunctions. This can 



be done at any T, at least numerically, using the implicit 
form of the exact fixed point given in the last section. At 
low temperature it can be done analytically, provided one 
distinguishes the two regimes, u ~ T (TBL) and u ~ 1. 
In the second regime the analysis becomes identical, to 
leading order in T, to the one performed directly at zero 
temperature in section III CI In the TBL, it is natural 
to look for solutions of the form f a (u) — f(u/T). The 
matching between the two regimes will be studied in the 
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next section. Inserting in the linearized version of eq. (|40[) 
and (|41[) and using the low T expansion of eq. ([55]) we 
obtain: 







(e 



- 2a) f{x) 
e 2 d 2 



(42) 



T0 2 (*))/>)] 



from which a systematic low T expansion of f(x) can be 
obtained. 

For a = e/3, one finds that a simple ansatz almost 
solves it (in an approximate sense given below) namely: 



1 



12T 



4>{x) 

noting that, from eq. ((38]) : 



ln(cf>(x)). 



4>'{x) = 
4>"{x) = 



X 

360 
1 

360 



£0' 



360 3 



(43) 

(44a) 
(44b) 



inserting into the right hand side of the above equation 
leads to: 



12T 

= (l-2/3)(l/0 ln< 

e 

it d 



(45) 



12 

30 



'(1 + In0) + ^- 



'ifl +0(T) 



= ^-^ + w +e {j) +0{T) 

Hence the dominant terms cancel, and the remaining 
term is found to be subdominant at large u going 

as ~ as u — > oo. A more complete analysis is done 
in section fV A[ where it is shown that the term e(02/0)" 
is also subdominant. For now suffice it to note that the 
form of eq. (I43[) is exact to dominant order in T for all 
x, and to next order in T it reproduces the exact large 
x behavior. Thus, surprisingly, we find evidence for a 
logarithmic correction to the eigenfunction, emerging at 
non-zero temperature T > 0. 



C. Logarithmic temperature dependence of the 
eigenvalue 

Let us examine the correction, induced by tempera- 
ture, to the eigenvalue. Assume that the main contri- 
bution to the eigenvalue correction comes, to first or- 
der, from the much larger regime outside the boundary 
layer where the eigenfunction is given by the expression 
of eq. (j2"3")l . Assume further that a is close to the value 
a = e/3 and expand the expression for the corresponding 



eigenfunction in powers of 5a 
u = 0, using eq. (j!M|) : 

1 



i, as well as around 



f a {u) = l + 5a ~-3-21n(4u 



+ 0(5a 2 ,u 2 ,5au) , u = O(l) 



(46) 



The divergence at u = is rounded off inside the TBL. 
The previous paragraph gives the expression inside the 
TBL, to leading order in T: 



/«(«) 



5a i 



1 



6T 



0(T) (47) 



c a \ 
6T J 



where we have multiplied with a constant in order to 
match the 1/u term of eq. (|4"6"|) for large u/T. 

The eigenfunction is now integrable and we enforce the 

1/2 

condition of zero mean, expressed as J Q duf a (u) = 
when taking into account that the eigenfunction is sym- 
metric about u — 1/2. We split the integral into two 
parts, inside u < kT and outside u > kT the boundary 
layer (any large constant k will do), and we find up to 
order 0(T, 5a): 



= 



duf a (u) 



(48) 



o 

5a e 



kT 



du 



du 1 



kT 







(—) 

3-21nAu + 0(u] 
5a (lnT + 0(l)) + 0(T). 



d„ [ - 

u 



This shows that for the eigenfunction to integrate to zero, 
i.e. for the InT divergence to be compensated, the eigen- 
value acquires a logarithmic temperature dependence. 



5a 



1 



21ni 



0{5a 2 ) 



(49) 



In fact, inserting this expression for the eigenvalue correc- 
tion back into the eigenfunction expansion of eq. (|46[) and 
normalizing by 6T/5ae, we retrieve exactly the asymp- 
totic form of the eigenfunction inside the boundary layer 
(eq.®). 



(50) 



6T , . . 6T 12T, u 

j—fa{u) = ln- 

5ae eu e 1 




12T 



m,/i + (- 



including all logarithmic terms. Hence our solution sat- 
isfies the required conditions, zero mean and matching 
between inside and outside the TBL. Note that it was 
necessary to not only include the dominant contribution 
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in the TBL, but also the subdominant one to logarithmic 
accuracy. 

In principle the low T expansion can be pursued to 
higher orders. The second order 0(Sa 2 ) is much harder 
to calculate though, because it requires us to exactly in- 
tegrate the non-divergent part of f a (u) over the whole 
interval [T, 1/2]. 



IV. NUMERICAL STUDY OF THE RANDOM 
PERIODIC CLASS AT T > (0 = 0, d = 2) 

To test the subtle mechanism for the selection of the 
eigenvalue based on the boundary layer matching, we now 
turn to a numerical analysis. 



A. Shooting to solve eigenvalue problem 

We numerically solve the FRG flow equation for the 
two-copy correlator (eq. ([ID])) linearized around the exact 
implicit solution of the one-copy correlator at finite T 
(eq. Ij31|l ). The eigenvalue equation is 







(e - 2a 



A* T "(u))f a (u) 



2A* T '(u)f a (u) 



(T + A* T (0)-A* T (u))f':(u) 



(51) 



with periodic boundary conditions: / a (0) = / a (l) = 1 
and /„(0) = f' a {X) = 0. Again, since the eigenfunction 
is even about u— 1/2, it suffices to consider the interval 
u E [0,1/2], requiring f^l/2) = 0. We set e to unity 
since it plays no role. 



1. Numerical details 

We need a continuous numerical representation for the 
fixed point correlator A^(u), whose analytical properties 
are given in section IIII Al The first step is to evalu- 
ate A^(0) for given temperature T, from the constraint 
1/2 = ^372 G(T, y + ). We have to find A^(0) such that 
the elliptic integral is exactly one half. This is easily done 
using the Brent method. The integrand has a divergence 
~ l/y/u at each limit, which is handled by a change of 
variables to v 2 — u. We use a Romberg integration rou- 
tine to efficiently get the desired precision. 

We then use the implicit expression of eq. ([3"T]) to cal- 
culate a discrete representation {itj(Ai), A,; i — 1..N} of 
the correlator. The Aj are chosen such that we have a 
sufficiently fine discretization inside the boundary layer 
u ~ T, even for very small values of T, and are chosen 
less dense outside the boundary layer, to reduce the total 
number N of support points. 

A cubic spline interpolation of the discrete function 
Aj(uj) allows us to obtain a continuous representation 
of the correlator Aj,{u). At the origin u = we take 
advantage of the exact series expansion of the correlator 



of eq. (|34]) . matching the series expansion to the spline 
at about u ~ 0.01T. 

With this continuous numerical representation of the 
correlator fixed point, we can solve the eigenvalue prob- 
lem of cq. (|5ip to arbitrary precision, using the shooting 
method: make an initial guess for the eigenvalue a, and 
integrate the ordinary differential equation starting from 
the initial condition / o (0) = l,/^(0) = 0, by means of 
a standard integration routine (e.g. odeint). Aiming for 
periodic boundary conditions f' a (l/2) = 0, one finds the 
eigenvalue a(T). 

In this way, we calculate the eigenvalues a(T) and 
eigenf unctions f a (u) as a function of the temperature T 
to arbitrary precision. We need however to go to very 
small values of T ~ 10~ 16 , to clearly see the logarith- 
mic dependence on T of the eigenvalue (eq. This 
requires quadruple precision for the numerics. 



B. Numerical results 

Indeed, we find numerically that the largest physical 
eigenvalue | is equal to one third plus corrections loga- 
rithmic in the cutoff T (see figure [TJ upper half). This 
confirms our analytical finding of logarithmic corrections 
to the eigenvalue. There exists one larger eigenvalue 
- = 1/2. It is easy to see that the exact eigenfunction 
for this eigenvalue is: 



fa- 



K 



»/2l 



A* T (0) - A* T {u) + T 



(52) 



for all u inside and outside the TBL (which correctly 
matches the eigenfunction f a =e/2 ~ l/( w (l — u )) f° r u ^> 
T). This eigenvalue does not acquire any corrections in 
the cutoff T, and most importantly, its corresponding 
eigenfunctions are strictly positive, they do not have any 
zeros. This means they cannot have zero mean and hence 
do not correspond to a correlator of the RP class. 



We recall the first analytical terms of the eigenfunction 
of eq. {SUJ with e = 1): 



lim 

u/T^oo 



1 



12T \n\ 1 



1- 
Mu) 



36T 2 



(53) 



36T2 

6T 

u 



- 12T In ■ 



T 



This expression compares rather well to the numerically 
calculated solution (figure [3]). Moreover, setting equa- 
tion (|53|) to zero gives us the first order term of the zero 
of the eigenfunction uq = 2 \ni/T ' ^ e P re f ac t° r \ is 
exactly the one found in the numerical data of the first 
zero of the eigenfunction (see figure [1] lower half). Thus 
the numerical results confirm our analytical analysis, pro- 
viding us with a coherent picture of the solution of the 
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a(T) 




0.04 
0.02 
0.0 



~i 1 r 

0.05 0.1 

l/ln( 1/T) 



FIG. 1: Slow logarithmic dependence on cutoff T. 
Eigenvalue a and b) ( x ) zero uo of the first eigenfunction ap- 
proach their zero temperature values a) 1/3 6)0 as l/ln[l/T] 
(we set e = I). The linear coefficient of the fitted numer- 
ical eigenvalue data coincides with the theoretical value of 
— 1/2. Equally, the linear coefficient of the fitted zero of the 
eigenfunction is identical to the theoretical value of 1/2. The 
ordinate scale is linear in l/ln[l/T]. As a guide, the corre- 
sponding value for T is given on the upper ordinate scale. 



linearized FRG flow equation at non-zero T. This is rep- 
resented schematically in Fig. 



As the cutoff T approaches zero, the eigenfunction 
shifts more and more weight into the ever smaller BL, 
in order to still fulfill the zero-mean constraint J f = 0. 
While this picture is satisfactory for any fixed T it is 
still not clear whether it could help to solve the prob- 
lem directly at T = 0. The question of how in the limit 
T — > this eigenfunction develops the non-analyticity 
necessary to balance the cusp in the one-copy correlator 
is especially subtle given that there are two regimes with 
different T-scaling properties. The boundary layer dis- 
appears as u B l ~ T, whereas the zero of the eigenvector 
approaches zero as uq ~ l/ln[l/T] (figure [TJ lower half). 



C. Consequence of logarithmic correction 

The unusual logarithmic correction to the eigenvalue 
— caused by a finite cutoff length — implies that it is 
very hard to calculate the latter by means of intuitive 
numerical approaches. If for example one simply tries to 
numerically integrate the FRG flow equation, one neces- 
sarily introduces a cutoff length (~ iV -1 , N being the 
number of discretization intervals in real space, or the 
highest frequency mode in Fourier space). This cutoff 
length has exactly the same effect as the finite tempera- 
ture cutoff, i.e. preventing access to smaller lengthscales. 
Thus even at zero real temperature, any finite numerical 







\ 

< >\ 


/ 


\u ~ 1/ log 1/T 


—A 



u 



FIG. 2: Different regimes of eigenfunction. The eigenfunc- 
tion has three different regimes, separated by two crossover 
values, represented schematically here. Near zero, the eigen- 
function presents a boundary layer of width u ~ T as found 
in Fig. [3] Outside the boundary layer, the eigenfunction 
falls off, and crosses zero at uq. This zero has a logarithmic 
dependence on temperature, introducing a second crossover, 
which scales differently from the width of the boundary layer 
uo ~ 1/ log 1 /T. For even larger values of u, the eigenfunction 
flattens but remains negative, such that the overall integral 
of the RP eigenfunction is strictly zero. 



1.0 



fr( u ) 



0.8 



0.6 



0.4 



0.2 




10 



20 25 30 35 



u/T 



FIG. 3: Scaling form of eigenf unctions. Inside a boundary 
layer of width u ~ T, the eigenfunctions present a functional 
form independent of cutoff T. Numerical data for various 
temperatures T = 10~ 4 ...10 -8 is indistinguishable from the 
theoretical expression fx (it) = 1/ yl-j- u 2 /36T 2 . 



cutoff introduces an immediate and non-negligible cor- 
rection to the eigenvalue of order lrLfiV] -1 . For example, 
a reasonably large N w 10 6 leads to a correction of order 
10% to the eigenvalue. 

A logarithmic correction to an eigenvalue, not unsimi- 
lar to the present situation, has been noticed in the con- 
text of a propagating wave frontsa. 
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V. GENERAL ANALYSIS AT T > AND 
EXTENSION TO OTHER CLASSES 

We turn from the simple case (d = 2, RP class) where 
the temperature does not flow, to the general problem at 
non-zero temperature. The two equations (JSJ) are studied 
with the temperature allowed to flow, i.e. 8 > 0. Despite 
Ti = Te~ el flowing to zero, it is still possible to use linear 
analysis, as we will show. We look again for a zero mean 
eigenfunction, placing us in the RP class for d > 2, and 
the RB class. 



A. TBL for one copy correlator 

We start by solving more accurately the first equation 
i]) and write: 

A(tt) = A* (it) , u = 0(l) (54a) 
A(tt) = A(O)-T i 1 ( xu /r/)-r ; 2 2 (x"/^)(54b) 
+0{T?) , u = 0(Tl) (54c) 
A(0) = A*(0)-T /7 i + O(T 2 ) (54d) 



with T t = Te~ m . One must have 0„(O) = 0, and (f>i{x) ~ 
\x\/& at large x to fit the cusp, hence the choice x = 
6|A*'(0 + )|. The zero temperature FP A*(u) can now be 
of RP type, as well as RB (in which case £ is non-zero 
and determined by the FP equation and the SR boundary 
condition). From the analysis performed in Appendix [Bl 
one finds: 



<f>i(x) = a/1 + x 2 /36- 1 = (f>(x) - 1 (55) 

and x 2 = e 2 X 2 = 36A'(0+) 2 = 36(e - 2C)A*(0) ( X = 1 
for the periodic FP) using 0"(O) = 1/36. From which we 
recover the zero temperature fixed point (eq !14p in the 
large argument limit [4>i{x — > oo) = ^ + j|y + . . . ]. To 
next order one finds: 



<^2 (x) 



X 2 H X ) 



[12(6-2e + 2C)(#r)-l) (56) 



+ x 2 (l- 11 ) + lx 2 <f>(x)(C-e)) 
6 

-3(4 + C-e)^ArcSinh(-))l 

6 

Note that <j) 2 {x) ~ (C - e)x 2 /(6x 2 ) + 0(x) at large x. 
This is compatible with A*"(0+) = (e - 0/3. Note also 
that: 



lim 



h (x) 
<j>(x) 



72 

-5-^ [ 156» + 7e - 22C 

X xr 



(57) 



3(e + 26»-3C)ln(-) ) + 0{x 



as promised in section llllBl hence validating the approx- 
imate solution given there. 



B. Equation for the two copy correlator 

Now we define the solution of the second equation in 
§ to be: 



D(u) = A(u) + f(u) 



(58) 



and we study the resulting equation for f(u) in a linear 
approximation: 



ft/(u) = (e-20/(u) + Cti/'(u) 

+ ^ T l + A(0)-A(u) )/(«)] 



(59) 



The only neglected term is — ^-£jjf(u) 2 on the r.h.s. One 
can write: 

dtf(u)= (e-2C) f(u) + (uf'(u) (60) 
+ &[{ T^ixu/Tt) + TfMxu/Tt) 
+ 0{T?) )/(«)] , M ~T ; 

$/(u) = (e-2C) /(u)+Ct*/'(u) (61) 
+ £[( A*(0)-A*(u) )/(«)] 
+ 0(T ; ) , tt~0(l) 

where we have explicitly separated the inside from the 
outside of the TBL. 



1. Eigenfunction inside the TBL 



Inside the TBL we look for a solution of the form: 

/(«)4/> = Wri) . u ~ T i ( 62 ) 

6T; 

As we will see below, the prefactor 1/TJ is crucial for 
obtaining a correct matching to the u = 0(1) regime. 
Equation (|6"Tj) gives: 



dj(x) = ( e -2(-6)f(x) + (<;-9)xf , (x) 

, X 2 ^ r 



( 0(ar)+r,02(a:)+O(lf 



fi dx 2 

The solution seems to admit the expansion 

1 



(63) 



/(*) = e 



2d/ 



<f>{x) 



+ T l yj(x) + 0(Tf) (64) 



where: 



= (e-2C-2a 



^ 4>(x) ^ ^ 4>{x) 2 



dx 2 



This yields: 

ip(x) 



<h(x) 
cf>{x) 

<h{x) . 



■ 4>(x)ij)(x) 



(65) 



i ( 36(e + 6 - 4C - 2a) 



0(x) 2 x 

, 6xArcSinh(f 
+ (2a-e + 3C)- 



(66) 
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with a undetermined. At large x one has: 



such that A a 



0. This is equivalent to the shoot- 



36 ,„ 

— 2a -e 

X 

36, 
+ — e + i 

X 2 



3C) In 



3 

4C- 2a 



:(C- 



(67) 
0(1/*) 



If one defines the "eigenfunction" f a (u) through f(u) — 
e 2al fa(u), it has the following behavior at large u/Tf. 



/«(«) 



X 

6T; 



6T; 36 ,„ 
— + — (2a- 

X u X 



(68) 



consistent, up to a normalization with the result (|50p for 
the RP class in d = 2, setting £ = 0, x — e > a = 1/3- 
This expression will be matched to the small it behav- 
ior of the eigenfunction f a (u) in the regime u = O(l), 
which contains a small l/u term. Note however that 

/»(«) ~ (x/(6T/))/Vl + X 2 w 2 /(36^ 2 ) in the TBL, hence 
/ a (0) has a different dependence in I that f a (u) for 
u = 0(1), i.e. it is a non-uniform eigenvector. 



2. eigenfunction outside the TBL 

For it of order unity one must study the linear differen- 
tial equation (|62p containing the zero temperature fixed 
point A* (it). 

Let us start with the RP class (for d > 2). Since 
the fixed point A* (it) does not change form (apart from 
the overall factor of e already taken into account) we 
expect the same behavior as for d = 2. From the 
ansatz of eq. (|62|) one sees that for £ = the behavior 
(eq. (|68p) identifies with eq. (|5D|) which matches the one 
for it = 0(1), (eq. (|4*6")l ). as discussed already in section 
IIII CI Hence for the RP class^, we find the same growth 
exponent a = e/3 for the function outside the TBL (for 
u = O(l)). But from the discussion of the end of the 
previous paragraph, the growth exponent of: 



e, 



2al 



(69) 



is determined by /(0) inside the TBL, hence a = a + |. 
This non-uniformity is the main difference with the case 
d = 2. 

For the RB class the equation for the eigenvector for 
it = 0(1) reads: 



= (e-2C-2o)/(u)+C«/ / («) 

+ A*(0)-A» )/(«)] 



(70) 



This equation involves the RB fixed point, which is non 
trivial— and was determined with high accuracy ir*i£ to- 
gether with the value for ( = 0.208298063. An analy- 
sis near it = shows that any solution of eq. ([70)) has 
the form of eq. (pM]) at small it, with B a — (e — 3C — 
2a)/A*'(0 + ). We will assume that one can then pro- 
ceed as for the RP case and look for the value a = orb 



ing problem of solving eq. l[70|) imposing that /(0) = 1 
and that /(it) decay at infinity. This fixes a unique and 
non-trivial value for a^B- We have solved this shoot- 
ing problem numerically using Mathematical, and found 
orb = 0.083346(6)e. The corresponding eigenvector is 
everywhere positive and integration of eq. (I70p easily 
leads to the following constraint: 

/>oo 

(e-3C-2o)/ duf(u) = -A'(0+)/(0). (71) 

Indeed one finds a^B < ( e ~ 3£)/2 (the + means that 
the domain excludes the TBL). Once this eigenfunction 
is determined in the region u = 0(1) the method to sat- 
isfy the zero integral condition over the full axis includ- 
ing the TBL is the same as for the RP class. First, one 
checks that the solution (eq. (|6"B"|) ) in the TBL matches 
correctly the small it behavior (eq. ([24]) ) of the it = 0(1) 
regime, using B a — (e — 3( — 2a)/A*'(0 + ). Second, one 
has again A aRB+ s a ~ Sa and one can proceed as in sec- 
tion [11121 One finds that the zero mean condition again 
leads to logarithmic corrections ai = a^B—K/ln(l/Ti) = 
orb ~ K/(0l) where K = duf(u) /d a A a \ a=ans given 
by eq. (f7Tj) . The main difference with the RP class, 
besides the value of the growth exponent, is that the 
logarithmic temperature corrections are also (weakly) l- 
dependent^. 



VI. SUMMARY AND DISCUSSION 

In this paper we have studied the problem of two mani- 
folds pinned in slightly different random potentials at the 
same temperature. We have written the coupled FRG 
equations for the single and two-copy correlators, A(it) 
and D(u), and temperature, to one loop accuracy. We 
have investigated how the difference between the two 
copies increase with scale. We have focused on zero- 
momentum (g = 0) quantities and specifically we have 
computed the fluctuations of the difference in (thermally 
averaged) center-of-mass positions U l (v), i = 1,2, of the 
two copies, in the presence of a uniform confining har- 

centered at 



-2l\ 



monic potential (of curvature m 
a common position u = v. This observable is exactly 
given by: 



(fii( w )-u2(0)-«) 2 (72) 
= A d (LTO)~ d TO~ 2C (A(i;TO C ) - D(vm c )) 



and can be seen to measure the r.m.s shift in position 
of (L/L m ) d roughly independent pieces of manifolds of 
typical size L m ~ \jm (hence the factor (Lm)~ d from 
the central limit theorem). The deviation of the one- 
copy center of mass from the center of the well is typically 
U % (v) — v ~ 0{m^^). We have defined the two growth 
exponents: 



0, := A(0)-D(0) 
A(u) - D(u)\ u=0{1) 



Ce 2hl 
Ce 2al 



fa{u) 



(73) 
(74) 
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allowing for the possibility that a ^ a. The coefficients C 
and C vanish when the difference in disorder 6 between 
the two copies is taken to zero. 

We have obtained these exponents at non-zero tem- 
perature, using linear analysis and matching the regime 
u = 0(1) to the regime u ~ T; = Te~ 01 in eq. ([74]) . 
corresponding respectively to v = 0{m~^) and v = 
0(Tm 9 - c ), the width of shocks, in eq. (|72|) . For the ran- 
dom periodic class for d > 2 and for the random bond 
class we found, with e = 4 — d: 



a+ 2 



a RB 



3 21n(l/T z ) 
0.083346(6) 



Ml/T,) 



(75) 
(76) 

(77) 



The result for the d — 2, RP class has been thoroughly 
checked via a numerical simulation. The other results 
assume that the scenario demonstrated in that case can 
be extended, which appears to be consistent. To confirm 
it further would require extensive numerics or a more 
complete analytical study. 

Since it originates from linear analysis, the growth 
of eq. (|73[) is valid only for a limited range of scales 
/ = ln(m/mo) (i.e. of masses m in eq. (|72p). From 
eq. (p)5|) , where the term (f 2 )" has been neglected the 
linear analysis is valid only, in the regime u = O(l), as 
long as Ce 2al <C 1. Beyond that scale the growth is 
non linear. However, since the eigenvector is peaked in 
the TBL the condition of validity of our analysis is more 
stringent. One can write the condition (f 2 )" <C Tif" 
and substitute / ~ (C/u)e 2al with u ~ T ; . This yields 
Ce 2al /Ti ~ Ce 2al -C Ti which, not surprisingly, can be 
written as: 



e, « t x 



(78) 



i.e. the width of the chaos BL (which exists at T = 0) is 
smaller than the TBL. This is the condition for validity 
of the thermal regime and linear analysis. Qualitatively, 
it means that the thermal width of the shocks should be 
larger than the typical shift in their relative position in 
the two copies. 

More issues remain to be understood. If one defines 
the overlap length L$ ~ e ls from 0; a ~ 1 one can ar- 
gue that the chaos exponent is a — a in the thermal 
regime. Note however that this assumes that C ~ <5 2 , 
a natural condition, but which may be spoiled at scales 
around the Larkin length where we do not have good 
control on the flow. Concerning the definition of Ls from 
eq. ([5]), using eq. (fTl"]) as was discussed there, we need a 
more precise calculation of the scale dependence of the 
non-local terms in the FRG. Finally, the present analysis 
does not solve the question of the zero temperature chaos 
boundary layer when the condition of eq. (|78p is violated, 
although it gives some insight into selection mechanisms 
for growth eigenvalues in such FRG equations. 



To conclude, we have made a step towards solving the 
intricate non-linear coupled FRG flow of the force cor- 
relators. We have found an interesting result for the 
growth exponent of the elastic manifold in frozen dis- 
order when the problem can be solved using a linearized 
flow equation. We have shown that it acquires a sur- 
prisingly large logarithmic correction when a finite cutoff 
length is present (here the temperature). More work is 
necessary to completely solve the problem and to estab- 
lish the exact relation between the eigenvalue found here 
and the chaos exponent. 



APPENDIX A: ANALYSIS OF LINEARIZED 
EQUATION 

We give a detailed account of the linearized flow equa- 
tion at zero temperature (eq. (|2"2"|) ). As we are looking 
for a continuous solution periodic on u € [0,1], we are 
only interested in the even solutions. The first solution 
that springs to mind is f(u) = l/u(l — u), corresponding 
to the eigenvalue a = e/2. It is however not integrable at 
u = 0. 

In order to get an idea of possible integrable eigenfunc- 
tions, we look for finite series solutions. Applying the 
Frobenius method, we postulate a power-series solution 
of the form 



/»(«) 



OO 

u ) j c n u 

n=0 



(Al) 



insert it into equation (|2"2"|) and equate coefficients of each 
term in the power series. The lowest power (indicial equa- 
tion) provides the two possible values for r = 0,-1. 

In the case of r — 0, the recurrence relation for the 
coefficients is: 



c m +i 



m(m + 3)-4 - 12a/e 
'' (m + l)(m + 2) 



(A2) 



Notice that the ratio of coefficients tends to one 
lim m _ >oo " +1 = 1, hence the infinite series does not con- 
verge uniformly on the interval [0..1]. We thus have to 
demand that the series terminate after a finite number 
of terms, i.e. that the m th coefficient be zero c m +i = 0, 
giving us the eigenvalues d m : 



1 



m(m + 3) 



(A3) 



Here we include even and odd solutions, whereas the 
eigenvalues a n of section III CI are a n = a m=2 n ■ The case 
r = — 1 does not lead to an independent solution, but 
points us into the right direction. Noting that the deriva- 
tive of In is given by — 1/u ( 1 — u) , a term that cancels 
the polynomial occurring under the second derivative of 
equation |22|) . we insert the following ansatz 



/-(«) 



q(u) 



u(l — u) 



p(u) In 



1 - it 



(A4) 
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with p(u),q(u) polynomials in u. Equating terms in 
In as well as terms in l/u(l — u) leads to the fol- 
lowing two equations: 

= 4(1 - 3a/e)p(u) + 2(1 - 2u)p' (u) (A5a) 

+ u(l -u)p"(u), 
= 6(1 - 2a/e)q{u) + u(l - u)q"(u) (A5b) 

- (1 - 2u)p(u) - 2u{l - u)p'{u). 

p{u) obeys the same equation as f a (u), hence the second 
type of logarithmic solutions has the same set of eigen- 
values as the first type (eq. (|A3|) ). and these eigenvalues 
(including a = e/2) are the only ones possible. 

p(u) gives us q(u), again by comparing coefficients, and 
thus we have the complete set of finite-series solutions. 
The first few are listed in table HI None of the finite-series 
eigenfunctions, however, fulfill the zero-mean condition 

duf a (u) = 0. From equation (|2"2")l one has 



APPENDIX B: THERMAL BOUNDARY LAYER 



duf a (u) = -—- — -j- . 

o 6 (! - 2a/e) 7o 



du[u(l - u)f a (u)]". 

(A6) 



In the case of the first type of solutions 



/ dtt/ (u) = \ - r [f a (0) + f a (l)], 

Jo 6(1- 2a/ e) 



(A7) 



which is non-zero for the even solutions we are interested 
in. 



In the case of the logarithmic solutions (eq. (|A4 



J duf a (u) 



6(1 - 2a/e) J 



q(u) (A8) 



1 — u . . 

+ u(l — u) log p(u) 

u 



6(l _ 1 2o/e ) [ q'(0)-q'(l)+p(l)-p(0)}, (A9) 



which is not zero either, for the even solutions. 



We now insert the ansatz of eq. (|54a[> into the FRG 
equations, for di(A(u) — A(0)) and SjA(O) respectively, 
collecting orders in 0(Ti), and using di = —OTiOt, one 
finds to first order: 



1 



= #'(0)- 
O=( e -2C)A*(O)- X 2 0'i'(O) = O 



(Bl) 
(B2) 



The solution is given in the text. The second order in T; 
yields: 

(0-Qx<fJ 1 (x)-(0 + e-2C)Mx)= (B3) 
-(0 + e -2C)7i=X 2 ^'(O) (B4) 



whose solution is: 



cj) 2 {x) 



4>{x) 
1 

+ — 

X 

1 2 

+ 2 X 



-4>'i (0)x 2 (B5) 

(12e-48C + 360)(0(x) - 1) 
6 + e-20 + ~x 2 ct>(x)((-e) 



3(e-3C + 20)a;ArcSinh(-) 



Fixing the free parameter (^'(O) by means of eq. (|B4[) . 
the above simplifies into: 



1 



x 2 4>{x) 



12(e-4( + 30)(<f>(x) - 1) (B6) 

' + e-2C)(l- 7 i) + ^ 2 ^)(C-e) 
6 



3(e-3C + 26»)a;ArcSinh(-) 



This can be further simplified using 8 = 2 — e + 2£ to yield 
the formula given in the text. One can check that this 
agrees with the result o£i*i once corrected for a misprint 
{f^>{ x ) ~> fs( x ) + x 2 (e — C)/3 in the formula given in 
Appendix D3 o&2). There the TBL expansion of the 
quantity y(u) = (A(«) - A(0) - Tf = T 2 f 2 ( X u/6) + 
X 2 T 3 f3(xu/6) + ■■ was computed. 



reven / \ 



/odd ( „ \ 
m m 



1 

2 
1 
:i 



_ J_ 

12 



Ce 



C ° (u(l 



Ce (l - 5u(l - u)) Co ( (1 ' 2 " ) l 1 a -°f~" )) ~ 12(1 ~ 5u(l - u)) In i=*) 



-u(l — 1() 

^ + 2 In i^i 

c D (l - 2m) 

)) 



TABLE I: First few eigenfunctions of the linearized flow equa- 
tion (1221) at zero temperature. 
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Strictly speaking, for the observable of eq. (J6]l one has the 
initial condition of the flow A(u) = A°(u) for m = +oo, 
i.e. I — — oo. 

As was shown there, such a term can be generated at T — 
only from the non-analyticity in the disorder correlator. 
Since the two copy disorder correlator R12 remains ana- 
lytic, such a term should not be generated in the two copy 
mutual correlation C12 which hence should remain a sin- 
gle logarithm. This is consistent with numerical results (G. 
Schehr, H. Rieger unpublished). 

This agrees with a droplet estimate in the case of two al- 
most degenerate wells in each copy, each thermally occu- 
pied with probabilities pt and 1 — Pi respectively, provided 
the probabilities of joint thermal occupations are products 
of the form pip2- It is obvious that in a given disorder 
sample (Vi, V2) thermal occupations are statistically inde- 
pendent in each copy if there are no interactions between 
the copies. However, there are obviously correlations be- 
tween the random variables pi and P2 with respect to the 
measure on the disorder. 
That is R'(0 ± ) = 0, i.e. no supercusp. 
Note again for £ = the exact eigenfunction of (|70p 
fa=e/2(u) = A * (0) f A . w which does not satisfy the re- 
quired conditions. 

We thank Kay Wiese for providing high precision approx- 
imations for A*. 

This dependence being much weaker than the leading ex- 
ponential one, it implies an additional I dependence of the 
eigenvector which is subdominant and neglected here. 



